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Given an assembly of viscoelastic spheres with certain material properties, we raise the question 
how the macroscopic properties of the assembly will change if all lengths of the system, i.e. radii, 
container size etc., are scaled by a constant. The result leads to a method to scale down experiments 
to lab-size. 
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Assume the dynamics of a certain granular system S 
is known. Will the dynamics change if we rescale all 
sizes by a constant factor a, i.e., R[ = aRi, but leaving 
the material properties unchanged? Q If scaling affects 
the system properties, how do we have to modify the 
material properties to assure that the system S and the 
scaled system S' behave identically? 

It is frequently desired to investigate large scale phe- 
nomena in granular systems experimentally, such as geo- 
physical effects or industrial applications. To this end 
one has to rescale all lengths of the system to meet the 
restrictions of the laboratory size, i.e. big boulders in the 
original system become centimeter sized particles in the 
experiment. Of course, one wishes that the effects which 
occur in the large system occur equivalently in the scaled 
system too. With the assumption of viscoelastic particle 
deformation we will show that naive scaling will modify 
the properties of a granular system such that the original 
system and the scaled system might reveal quite different 
dynamic properties. To guarantee equivalent dynamical 
properties of the original and the scaled systems we have 
to modify the material properties in accordance with the 
scaling factor and we have to redefine the unit of time. 

As another consequence of the scaling properties we 
claim that for numerical simulations of granular material 
it is not sufficient to provide relative data such as, e.g. 
to describe the container size in units of the particle di- 
ameter. We will show an example where the dynamics of 
a granular system changes significantly with system size, 
although all relative sizes are kept constant. 

In a simple approximation, a granular system may 
be described as an assembly of spheres of radii Ri, 
i = 1, . . . , N. If two particles of radii Ri and Rj at posi- 
tions fj and fj touch, i.e., if £y = R4 + Rj — |fj — fj\ > 0, 
they feel an interaction force 



The normal force F n can be subdivided into elastic 
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and dissipative parts F n = F^\ + F c " is 



The elastic 



force for colliding spheres is given by Hertz's law 
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with R cS = R i R j /(R i + R 3 ) and Y, v being the Young 
modulus and the Poisson ratio. Equation ^ also defines 
the prefactor k which we will need below. 

The formulation of the dissipative part of the force F£- 
depends on the mechanism of damping. Here we will fo- 
cus on viscoelastic damping which is the most simple as- 
sumption for dissipatively colliding bodies Q. It implies 
that the elastic part of the stress tensor is a linear func- 
tion of the deformation tensor and the dissipative part 
of the stress tensor is a linear function of the deforma- 
tion rate tensor. It is valid if the characteristic velocity 
(the impact rate g) is much smaller than the speed of 
sound c in the material and the viscous relaxation time 
T v i s is much smaller than the duration of the collisions 
t c 0. The range of the viscoelastic model is, hence, lim- 
ited from both sides: the collisions should not be too fast 
to assure g <C c, t v - 1s -C t c , and not too slow to avoid 
influences of surface effects as adhesion. For viscoelasti- 
cally colliding spheres the dissipative part of the normal 
force reads ||, |7| 



A dt de el 



3 a r? d t 
= 2 Ak ^ Tt 



(3) 



The dissipative material constant A is a function of the 
viscous constants ?/i/ 2 , the Young modulus Y and the 
Poisson ratio v (for details see Q): 



A=- 
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with 



the unit vector in normal direction 
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and the respective unit vector in tan- 



gential direction ty. Eventually, external forces as, e.g. 
gravity, may also act on the particles. 



Combining the forces 
tion of motion 
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and (p|) one obtains the equa- 
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with m cff = mirrij / (rrii + Trij) and with the initial condi- 
tions £| t=0 = and d£/ dt\ t=0 = g. 
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In dimensionless variables, £ = £/£o, t = t/ro, Eq. (||) 
reads [pi 



dr 2 4 s 2 V 4 
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with £ =0 and g?£ / dr = 1 . The characteris- 

tic length £o 1S the maximal compression for the equiva- 
lent undamped (elastic) problem which can be found by 
equating the kinetic energy of the impact m g 2 /2 and 
the elastic energy at the instant of maximal compression 
2k£q^ 2 /5. As characteristic time To we define the time in 
which the distance between the particles changes by the 
characteristic length just before the collision starts 
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t = Zo/g. 
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The only term in Eq. (|6j) which depends explicitly on the 
system size and on material properties is the prefactor 
in front of the third term. If the scaling procedure af- 
fects the value of this term it will change the dynamics 
of the system. To ensure identical behavior of the scaled 
system, however, besides the identity of this prefactor, 
further requirements have to be met which will be dis- 
cussed below. 

Expanding our abbreviations we obtain 
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with p being the material density. The function 
F(Ri,Rj) collects all terms containing Ri and Rj\ 



F(Ri,Rj) 



[RiRj/ (Ri + Rj)] 
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Scaling the radii by a, the function F scales F(R' i ,R'A = 



F(aR iy aRj) 



' 1 F(R i ,Rj). Obviously, simple scal- 



ing of the system in general affects the prefactor Eq. (g) 
already via the scaling properties of F(Ri,Rj), hence, 
in general the original system and the system where the 
lengths have been scaled by a factor a differ in their 
dynamic properties. More explicitly, one can show that 
naively scaling the system by a factor a < 1 will lead to 
a comparatively more damped dynamics. 

To provide equivalent dynamical properties of the sys- 
tems, therefore, we have to modify the material proper- 
ties in a way to assure that the equations of motion of 
both systems are equivalent, which in turn assures that 
the prefactors (||) of the original system and the scaled 
system are identical. 

One of the few things which cannot be modified in ex- 
periments with reasonable effort is the constant of gravity 



G. That implies that going from S to S' not only G but 
all other accelerations must remain unaffected too 
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yielding t' = yfat. Hence, scaling all lengths x' = ax 
implies that times scale as t' — \foct if we require that 
the gravity constant stays unaffected. Thus, the clock in 
the scaled system S' must run by a factor ^fa faster (or 
slower) than the clock in the original system. In other 
words, if in the original system we observe a phenomenon 
at time t = 100 s, we will find the same effect in the 
scaled system at time t' = yfa 100 s. Scaling of time is a 
direct consequence of scaling the lengths if the constant 
of gravity has the same value in both systems. 

In the scaled system the equation of motion of a par- 
ticle contact reads 

d 2 e , «' ,,,,3/2 L 3 . k' m dl? _ 
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If we apply our scaling relations which were introduced 
above, i.e. 



we obtain 
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Comparing with Eq. (||) we find the conditions to assure 
identity of the equations of motion: 
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Using the definitions of k (Eq. (^)) and m cff yields finally 
Y V Y 
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A' = Va~A. (15) 



If we choose material constants which obey Eqs. ( |15| ) we 
will obtain the original equation of motion after scaling 
the system back to its original size, i.e., both systems are 
equivalent. 

When we incorporate the tangential force F* into the 
analysis, of course, we have to require that this force 
scales in the same way as any other force, namely 
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given that accelerations are invariant under scaling. This 
requirement has to be met by appropriately scaling the 
material constants, particularly the friction constant, re- 
sulting in an additional scaling equation. Its form de- 
pends on the underlying friction model, i.e. on the func- 
tional dependence of the tangential force on the geometry 



3 



TABLE I: Material parameters used in the simulations 



samples 


original 


scaled I 


scaled II 


plane length 


10 m 


62.5 cm 


160 m 


particle diameter 


10 cm 


6.25 mm 


160 cm 


P 


2 g/cm 3 


2 g/cm 3 


2 g/cm 3 


Y 


8 GPa 


0.5 GPa 


128 GPa 


A 


2-10~ 5 s 


5 • 10~ 6 s 


8- 10~ 5 s 


/' 


0.1 


0.1 


0.1 



and the material properties as well as on the compression 
and the relative velocity. For instance, if we assume the 
most simple tangential force law 

F* = fiF n , (17) 

with /i being the friction coefficient we can conclude that 
the friction coefficient has to be invariant with respect to 
scaling 

// = //. (18) 

In literature, one can find a variety of different models 
which might be more realistic for the description of the 
dynamics of a granular material (for an overview see, e.g., 
pot). Applying the same procedure for these laws will 
result in different scaling relations and, hence, in different 
conditions of the type Eq. (jT||) . A thorough discussion of 
the scaling of more realistic tangential friction laws will 
be published elsewhere. 

To demonstrate the derived scaling laws, we present 
results of a 2D Molecular Dynamics simulations of sta- 
tionary flow of 1000 particles down an inclined plane. 
In the "original system" the length of the plane is 10 m 
(periodic b.c), its slope is 15 degrees, and the particle 
diameters have been randomly chosen with an average 
of 10 cm. The material parameters are given in Table Q. 
Particles are subjected to normal viscoelastic forces as 
well as tangential forces (of the simple type discussed 
above) which, although not very realistic, is enough for 
the purpose of showing the scaling properties. Figure |l| 
shows a snapshot of the flow. 




FIG. 1: Snapshot of the granular flow down an inclined plane. 

The solid curve in Fig. ^ represents the density pro- 
file perpendicular to the inclined plane of the original 



system and of the correctly scaled systems I and II, tak- 
ing into account the necessary modification of the mate- 
rial parameters according to Eqs. ( p^| , ^8| ). The modified 
material parameters are summarized in Table |. It is 
not surprising that all three lines collapse perfectly, since 
in scaled variables the equations of motion of the origi- 
nal system and the correctly scaled system are identical. 
Therefore, particle trajectories and, hence, the density 
plot must be identical. 
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FIG. 2: Density profile of the flow down an inclined plane. 
Thick solid line: original system together with scaled systems 
I and II (ai = 1/16, an — 16, correct scaling, parameters 
see Table |j) . The curves collapse precisely to one single line. 
Dashed line: scaled system II, a = 1/16 (incorrect scaling), 
dot-dashed line: scaled system, a = 16 (incorrect scaling). 
Thin solid line: system incorrectly scaled by a — 49. For the 
scaled systems, the horizontal axis has been multiplied by 16 
(I), 1/16 (II) and 1/49. 

Figure |^ also shows two curves for the stationary den- 
sity profile of the systems scaled by a = 1/16 (system I) 
and a = 16 (system II) where only the lengths have been 
scaled but erroneously the material parameters have been 
kept constant (incorrect scaling). These curves deviate 
from the profile of the original systems since the original 
and the scaled systems obey different equations of mo- 
tion. As anticipated, the scaling by a — 16 only applied 
to lengths leads to a more elastic dynamics, i.e., com- 
paratively more dilute flow (dot-dashed curve in Fig. ||) . 
This effect is amplified when larger scaling factors are 
applied (a = 49), where in twice the simulation time, 
no steady state is found while the density profile keeps 
relaxing towards more dilute states. The instantaneous 
density profile is shown as a thin-solid curve in Fig. |^, as 
an indication that incorrect scaling can ultimately lead 
to completely deviated results. 

The simulations show that the simple length scaling of 
a granular system by a constant factor a changes the dy- 
namical properties significantly, if the material parame- 
ters are kept constant. In order to obtain identical results 
for a scaled system one has also to modify the material 
constants and redefine the unit of time. These necessary 
scaling relations are summarized in Table [Tl| 

The knowledge about these scaling relations offers the 
possibility to scale down real world systems, e.g. geo- 
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TABLE II: The necessary scaling relations when transiting 
from a granular system (S) to a scaled system (S") having 
identical dynamic properties 





original system scaled system 


all lengths 
time 

elastic constant 
dissipative constant 


x ax 
t y/at 

p(l-S) %(W) 
A ^a~A 



physical or industrial granular systems, to sizes where 
laboratory experiments can be performed. If one scales 
down such a granular system one has to replace the orig- 
inal material by a material which meets the scaling re- 
quirements discussed in the text. 

We want to give an example: Assume in the origi- 
nal system one deals with steel spheres (Y — 20.6 • f0 10 
Nm -2 , v = 0.29 and p = 7,700kg m -3 ) of average ra- 
dius R = 10 cm and system size of L = 10 m. The prop- 
erty whose scaling behavior is known is Y/{p{l — v 2 )) = 
2.92 • 10 7 m 2 s~ 2 . One wishes to know (to measure) a 
certain value at time t = 100 s. In the lab we per- 
form the experiment with an equivalent system of size 
L' = 1 m, i.e. we scale the system by the factor a = 0.1, 
including all radii. From the scaling relations we see 



that we have to find a material whose scaled property 
is Y'/(p'(l - v 12 )) w 0.3 • 10 7 m 2 s~ 2 . From tables |j we 
find that we can use plexiglass (Y = 0.32 • 10 10 Nm -2 , 
v = 0.35 and p = 1, 200 kg m~ 3 ) in order to obtain this 
value. Therefore, we have to perform the experiment 
with plexiglass spheres and have to measure the value of 
interest at time t' = 31.6 s. 

One can imagine that not for all scaling factors a one 
will find a proper material, however, nowadays it is possi- 
ble to manufacture materials which can meet demanding 
requirements, such as high softness along with a custom- 
designed damping constant. 

The scaling properties have also consequences for 
Molecular Dynamics simulations of granular systems: 
namely, it is not sufficient in simulations to provide rela- 
tive parameters such as the container size in units of the 
particle radius. As demonstrated, the result of a simu- 
lation (and, of course, also of a real world experiment) 
depends on absolute values. 
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